Sensitivity to the initial state of interacting ultracold bosons in disordered lattices 
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We study the dynamics of a nonlinear one-dimensional disordered system obtained by coupling 
the Anderson model with the Gross-Pitaevskii equation. An analytical model provides us with a 
single quantity globally characterizing the localization of the system. This quantity obeys a scaling 
law with respect to the width of the initial state, which can be used to characterize the dynamics 
independently of the initial state. 
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I. INTRODUCTION 



The Anderson model, introduced about 50 years ago 
is the simplest model describing the effects of disor- 
der in a quantum system in a relatively realistic way. It 
has however been recognized rather early that the model 
probably relied on too strong simplifications to match 
the actual behavior of electrons in a crystal: The An- 
derson model is a one-particle model at zero tempera- 
ture. Moreover, the difficulty of experimental investiga- 
tions in condensed matter systems prompted a search 
for other systems where such effects could be observed in 
more favorable conditions, leadingto studies of the local- 
ization of electromagnetic waves [1, 01 and sound waves 
0,f^, where particle interactions are obviously absent. A 
breakthrough has been recently realized by using laser- 
cooled atoms 04121 1 ^ development inscribed in a more 
general trend of studying man y- b ody systems with ultra- 
cold atoms in optical lattices [13fl . 

One of the long-standing questions concerning the An- 
derson model is the effect of particle-particle interactions. 
It was conjectured that the electron-electron coulomb re- 
pulsion would suppress localization. Including such ef- 
fects in the model, however, implies going from a "sim- 
ple" one-particle picture to a very complex many-body 
description. The problem of interactions is very rich 
in disordered ultracold atomic systems: For such sys- 
tems, the interaction strength can be driven from repul- 
sive to attractive, via the so-called Feshbach resonances 
[l^ . Moreover, by choosing the atomic species, the quan- 
tum statistics can also be changed from Fermi-Dirac to 
Bose-Einstein (see, e.g., [3| and references therein). One 
can conjecture, for example, that attractive interactions 
shall increase localization in a cold Bose gas, but not 
in a Fermi gas. All possible combinations of the sign of 
interactions and quantum statistics can be similarly con- 
sidered Mean- field theories, which simplify consider- 
ably the many-body problem, were generally considered 
to give a rather poor description of the electronic many- 
body problem in a crystal |17| . However, it became clear 
in recent years that these theories give, on the contrary, a 
very satisfactory description of ultracold Bose gases with 



(weak) interactions in a variety of situations \}^. This 
puts into evidence the interest of studying cold atoms in 
the presence of disorder. 

The first experimental result in this new field was the 
observation of "dynamical localization", a manifestation 
of the Anderson localization in the momentum space, by 
Raizen and co-workers, in 1994 [l^. Anderson localiza- 
tion of bosons in ID @ and ID localization in the Aubry- 
Andre model have also been observed experimen- 
tally, and^recently, observations of the 3D localization of 
fermions [l^ and bosons [ll| were claimed. The effect of 
interactions is also being studied experimentally (l9l-l2l| . 
Finally, the Anderson metal-insulator transition is being 
actively studied with a cold atom "quantum simulator" of 
the 3D Anderson model [22| , includingjthe experimental 
determination of its critical exponent [3, [23l . [2^ and the 
study of its critical state (25j . 



II. FRAME AND SCOPE 

We start by considering the tight-binding description 
of an one-dimensional (ordered) lattice of period d ob- 
tained by projecting the eigenstates of the spatially peri- 
odic Hamiltonian on a basis of localized functions, usually 
Wannier functions Wnix), associated with each site n: 



(1) 



which produces a discretized eigenvalue problem 



We use here the (usual) symmetric first-neighbors ap- 
proximation Tr ~ ~TSr,±i and, as usual, rescale time 
such that t — !■ Ht/T (and energies correspondingly: E = 
s/T) ill order to obtain the tight-binding equation: 



'^n^n ^n—l -^C^j 



(2) 



If all sites have the same diagonal energy, one can redefine 
the energy origin so that i)„ = 0, then (26l. [27|: 
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where q is the quasimomentum, E{q) are the energies on 
the first band, p is the density of states. The eigenfunc- 
tions tpq (x) are delocalized Bloch waves of quasimomen- 
tum q G [0, tt]. 

Tlie Anderson model postulates that the main effect 
of crystalline disorder is to randomize the Vn in eq. ([2|) 
in an interval [~W/2, W/2], an approximation known as 
"diagonal disorder". The introduction of disorder redis- 
tributes the eigenenergies between the values —2 — W/2 
and 2 -I- 14^/2 and, more importantly, produces a localiza- 
tion of the corresponding eigenfunctions (4^. We shall 
index the eigenfunctions byz^ = g — 7r/2so that = is 
at the center q = n/2 oi the energy band. The shape of a 
given eigenfunction depends on Ei, and on the realization 
of the disorder, but if one averages the eigenfunc- 
tions in a given energy interval [E^E + /S.E] over many 
realizations of the disorder one finds an exponential shape 

V..H-exp(^-^-^j 

(we use overbars to indicate averages over realizations of 
the disorder) where i^{W) is the localization length. In 
the ID case and in the limit of weak disorder one can 
show that (ill 

UW)^^cos':., (5) 

which vanishes at the borders of the band and has its 
maximum value £q ^ 96W~^ at the center of the band. 
This localization can be interpreted by noting that, in 
the absence of disorder, all sites have identical energies, 
so the particle can tunnel from a site to the next one, gen- 
erating a diffusive motion which leads, asymptotically, to 
a complete delocalization in the lattice. In the presence 
of diagonal disorder, however, it is unlikely that neighbor 
sites have close enough energies, so that in general, the 
particle can only perform virtual transition to neighbor 
sites, which leads to a localized exponential spreading of 
the eigenfunction. 

We take into account interactions by using a mean- 
field approximation. The decomposition eq. ([T]) is used 
to transform [ij the Gross-Pitaevskii equation (or non- 
linear Schrodinger equation) in a discrete nonlinear set 
of equations 

iCn = VnCn - C„_i - C„+i -|- g |c„|^ C„ (6) 

where g is the parameter characterizing atom-atom in- 
teractions, proportional to the s-wave diflFusion length 
(28I. [29j. Disorder is introduced by randomizing w„, as 
before. This model is occasionally named DANSE (Dis- 
crete Anderson Nonlinear Schrodinger Equation). As we 
use mean-field theory throughout this work, we shall use 
the terms "interaction" and "nonlinearity" interchange- 
ably. 

The main question in the problem of interacting dis- 
ordered systems is, will interactions destroy localization? 
The DANSE model has been used, e.g. by Pikovskii and 



Shepelyansky [sO], to investigate this question. They ob- 
served numerically a revival of diffusion, leading, at very 
long times (^ 10^), to a ^a;^) ^ sub-diffusive behav- 
ior with a Ri 0.3. Other works also tackled this prob- 
leni, with somewhat contradictory results. For example, 
in [3l|, the use of the so-called "participation number" 

1*^"!^) (instead of (x^)) as the quantity charac- 
terizing diffusion, lead to the conclusion that localization 
survives even in presence of a moderate nonlinearity, but 
subsequent studies j32l . |33| confirmed the existence of a 
subdiffusive dynamics, with an exponent a that typically 
depends on the initial state. Experimental evidence of 
subdiffusion in the Aubry- Andre model has also been ob- 
served recently [2l| . 

One can formally write the "diagona"l part of eq. ^ as 

(^Vn + <7|cn|^^ c„. which allows us to interpret the non- 
linear term v^^ = g|c„|^ as a "correction" (depending 
on the site population) to the energy of site n ^5|. This 
heuristic picture, although not rigorous, is often useful in 
interpreting the behavior of the system. 

Mathematically, much is known about eq. © in the 
absence of disorder, in different contexts (see e.g. 1341 ): 
As a special case of the Ginsburg-Landau equation [35| it 
can describe some aspects of superfiuidity, and, in optics, 
it describes the Kerr effect in a multimode system (36| . 
Two nonlinear effects play a particularly important role 
in our study. The first one is the called self-trapping, 
which manifests itself when a given site n has a much 
larger population than its neighbors. In such case, the 
correction v^^ decouples it from its neighbors (much as 
the disorder itself does), thus inhibiting diffusion. The 
second effect is the existence of chaotic behaviors (in a 
classical sense - that is, chaos related to sensitivity to ini- 
tial conditions) , due to the presence of nonlinearity (stI - 
I39I . The chaotic evolution of the amplitudes c„ may gen- 
erate strong variations of the nonlinear correction v^^ , 
eventually bringing, even if only for a short time, neigh- 
bor sites close to degeneracy, thus favoring diffusion. In 
the presence of both disorder and interactions, different 
regimes are possible: Localization is expected to survive 
(or at least be destroyed only at very long times) if inter- 
actions are weak enough compared to disorder, diffusion 
(or subdiffusion) induced by the chaotic evolution is ex- 
pected if these two effects are of the same order of magni- 
tude, and self-trapping is expected to dominate, inhibit- 
ing diffusion, in the regime of strong interactions. These 
regimes have indeed been observed numerically [s^, [13 ■ 
It can be tempting to classify these regimes by comparing 
to W, but one must not forget that, as depends 
on Cn{t), the dynamics depends on the particular trajec- 
tory of the system in the c„ space. Such a classification 
scheme is thus useless to define "phases" of the system. 

This is an important and - from the point of view of 
the familiar linear quantum systems - unusual character- 
istic of nonlinear systems: The initial state plays a very 
important role (much more than in linear systems) in de- 
termining the dynamics. For example, if one takes, as in 
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Figure 1: Survival probability p{g,t = 10^) as a function of the nonlinear parameter g for different values of the width Lo of 
the initial state and disorder (a) W = 2, (b) W = 3, {c) W = 4. Three regimes can be identified: Quasi-localized for low g; 
chaotic evolution with destruction of the localization for intermediate values of g; and self-trapping for high g. The nonlinearity 
destroys the localization almost completely in the chaotic regime, and self-trapping is more efficient, as expected, for smaller 
values of the initial state width Lq . For low disorder and small initial state width the localization due to self-trapping becomes 
much more efficient than the Anderson localization. Results are averaged typically over 1000 realizations of the disorder and 
of the initial phase distribution. Values of Lq: 3 (blue squares), 7 (green triangles), 13 (red diamonds), 21 (cyan stars), 31 
(magenta circles), 41 (yellow inverted triangles). 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0( 

1.0 
0.8 
0.6 
0.4 
0.2 
0.0, 





(a) 


♦ V 











5.10" 
t 



10^ 



(b) 



5.10" 
t 



10^ 



Figure 2: Comparison between numerical simulations (sym- 
bols) of p{g,t) and the fit (solid lines) by eq. ([8]) for Lo — 31 
and (a) g ~ 1, (b) g — 100. Values of W: 1 (blue squares), 2 
(green triangles) and 3 (red diamonds). 



[201, an initial condition c„(t = 0) = Sno, one strongly fa- 
vors self-trapping, and thus reduced diffusion. A choice 
like c„ — 1/Vn for |n| < N/2 favors, for moderate val- 
ues of N, chaotic dynamics, and thus diffusion. The main 
goal of the present work is to study the impact of initial 
conditions in the dynamics of the DANSE model. We 
shall show that this generally forbids a general classifi- 
cation of the dynamics in the parameter plane g, W . In 



the frame of a simplified model, introduced in sec. IIIIl 
we shall however show that one can define an "effective 
wavepacket length" obeying scaling laws with respect to 
the extent of the initial state (sec. lIVp. This allows us to 
define a generalized nonlinearity parameter, depending 
on the initial state, which leads to a more satisfactory 
classification of the dynamical behavior in the weak dis- 
order limit. Finally, we discuss in sec. |V] the impact of 
some important neglected effects on our results. 



III. THE MODEL 

The search for asymptotic behaviors of the type (a:^) ^ 
t" ill the DANSE model implies describing the system 
by an array whose size increases with the evolution, the 
number of elements typically increasing as This 
is not our aim here: If one starts with an initial state 
localized in a relatively small region, the restoration of 
(sub-)diffusive dynamics by the nonlinearity should be 
apparent on the fact that parts of the wavepacket must 
continuously "escape" the initial region. We thus study 
a disordered lattice in a box L containing L ~ 100 sites, 
starting from an initial state of width Lq. At the frontiers 
X — ±L/2 of this box we place "absorbers" by adding 
an imaginary part —irja to the potential, which increases 
exponentially over a distance La ^ 10, in order to prevent 
reflection of the parts of the wavepacket approaching the 
limits of the box. 

We are interested here in the effect of the spatial ex- 
tension of the initial wavepacket. In order to simplify the 
problem, we use a initial wavepacket of square shape and 
width Lq: 




-1/2 



|n| < (^0 - 
otherwise. 



l)/2 



(7) 



We prevent dominant quantum interference effects by set- 
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Figure 3: Effective localization length as a function of the nonlinearity g, for (a) W = 2, (b) W — 3 and (c) W = 4 and various 
values of the width Lo of the initial state. Note that for g ^ we retrieve the linear localization length £o ~ 96W~^ in the 
weak disorder limit W < 3. Same graphic conventions as in fig. [T] 
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Figure 4: Activation time t2 as a function of the nonlinearity 
g for (a) W = 2 and (b) W = 4 and various values of the 
width I/O of the initial state. Same graphic conventions as 

fig-ID 

ting random phases to the c„. We shall consider in seclVl 
the effect of non-random quantum phases. 

The part of the wavepacket that, at time t, has re- 
mained in the box L (that we shall call the "survival 
probability") is given by 

(L-l)/2 
n=-(L-l)/2 

We display in fig. mp{g,t) calculated a.t t ~ 10^ as a 
function of the nonlinearity g for various sizes of the ini- 
tial state Lq, for a disorder amplitude W ~ 2 (part (a) 



of fig. [T]), = 3 (b) and = 4 (c). The numerical 
integration is done by the Crank- Nicholson method [i^ ■ 
The probability density |c„p is averaged over typically 
1000 realizations of disorder {vn} and of the initial phases 
arg[c„(0)]. Note that one has p{g Q,t = 10^) < 1, due 
to the fact that the initial state typically projects onto an 
ensemble of localized eigenstates with localization lengths 
given by eq. ([5]) . Some of these eigenstates have localiza- 
tion lengths large enough to "touch" the borders of the 
box L, so that the projection itself leads to losses. One 
can see by comparing the three plots of fig. [T]that these 
losses decrease with the disorder W, and are independent 
of Lq, as expected for the linear regime. 

From fig. [T] one can identify three main dynamical 
regimes for increasing nonlinearity parameter g: i) A 
"quasi-localized" dynamical regime in which localization 
is only weakly perturbed by the nonlinearity; ii) a dy- 
namical regime in which localization is destroyed by the 
chaotic dynamics (iHl, the diffusion is reestablished and 
the losses increase notably; iii) self-trapping regime in 
which diffusion (thus losses) is inhibited again. It is inter- 
esting to note that, for concentrated initial states Lq < 20 
and weak disorder, self-trapping can be much more effi- 
cient than localization in suppressing diffusion. These 
dynamical regimes correspond to those studied e.g. in 
Laptyeva et ai. [ssj by considering the long-term behav- 
ior of {x'^y Fig.[T]shows that the crossover between these 
regimes strongly depends on the width of the initial state. 

In the absence of nonlinearity, the initial state evolves 
until it takes an exponentially localized shape of width 
~ £o and then "freezes". In the presence of the nonlinear- 
ity, one observes a global enlarging of the wavefunction, 
which results in a increasing absorption at the borders 
of the box L. In the weak disorder limit we can con- 
struct an analytic model for these losses, starting from 
the properties of the linear [g ~ 0) system. This model, 
discussed in appendix |X1 produces an analytical expres- 
sion describing the asymptotic behavior of the survival 
probability p{g = 0,t), eq. (|A5p . Remarkably, a small 
modification of this expression furnishes also a function 
describing the asymptotic behavior in the nonlinear case, 
for values of g as large as 10'^. This expression is 
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Figure 5: Scaling law. The same data of fig. [3] is plotted in terms of the scaled quantities g = gL'^^^'^ and I = £ct[LQ^^* for 
(a) W — 2, (b) W — 3 and (c) W = 4. A clear grouping of the curves in the nonlinear region is observed. Same graphic 
conventions as fig. [T] 
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Figure 6: False colors plot of the scaled effective wavepacket length ^ as a function of the scaled nonlinearity g and of the 
disorder W for (a) Lo — 21 and (b) Lo =41. The structure of the parameter space is very similar for the two values of the 
initial state width. 
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where ia is given by eq. (jA6p . depending on the geome- 
try and on the localization length (.^{W), and icS and t2 
fitting parameters. With respect to eq. (jA5[) of app. |^ 
we simply replaced £o by an "efi'ective wavepacket length" 
led which is the second fitting parameter of our model. 
Fig. [2] shows that the asymptotic behavior of the sur- 
vival probability is very well fitted by this formula. The 
quantity ^off plays a major role in the present work, as it 
can globally characterize the wavepacket shape. We at- 
tribute this somewhat astonishing property of our model 
to the fact that, even in presence of the nonlinearity, the 
wavepacket displays exponential wings in the asymptotic 
regime, as in can be verified numerically. 

We display in fig. |3] the behavior of the effective 
wavepacket length ^ett as a function of the nonlinear- 



ity, for disorder parameters W = \ [plot (a)], W = 2 
[plot (b)[, and W — Z [plot (c)[ and for various val- 
ues of the initial state width Lq. One clearly identi- 
fies the regions corresponding to the three dynamical 
regimes discussed above: Quasi-localized regime, char- 
acterized by a constant value of ^cff, independent of g 
and Lo! chaotic regime, characterized by a marked in- 
creasing of £eff with g and strongly dependent on Lq; 
and self-trapping regime, in which £off decreases again 
and can even become smaller than the Anderson local- 
ization length and which is also dependent on Lq- 
Fig. 0] shows the behavior of the "activation time" t2 (see 
App. |A| for the same set of states and parameters as 
fig. [31 This time roughly corresponds the time for the 
wavepacket to "touch" the border of the box L, when 
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losses become important. In the quasilocalized regime 
t2 is thus very high. The onset of the chaotic behav- 
ior, favoring diffusion, produces a dramatic decrease of 
t2 which becomes virtuahy zero (that is, losses began 
almost instantaneously). Self-trapping, inhibiting diffu- 
sion, produces an increase of ^2- The same three dynam- 
ical regimes can hence be observed also in the behavior 
of ^2- This fit parameter, although necessary to obtain 
a good agreement with numerical simulation, thus car- 
ries essentially the same information as i^s- For g ^ oo, 
one expects the wave packet to evolve very little due to 
"immediate" self-trapping, that is, its asymptotic width 
should be proportional Lq. In fig. [U £cS clearly has not 
attained this large g regime, but we observed numerically 
that even for moderate values of g ^ 300 it is roughly pro- 
portional to Lq: £cs{9 = 300) « 0.25Xo for W ranging 
from 2 to 4. 




Figure 7: Survival probability at t = 10^ for a square initial 
state (full lines) and a gaussian initial state eq. (|10p (sym- 
bols), for W = 1, Lo = 21 (blue/squares), W = 1, Lo = 41 
(green/triangles), W = 3, Lq = 21 (red/diamonds) and 
W = 3, Lo = 41 (cyan/circles). The dependence on the shape 
of the initial state is very small. 



IV. SCALING LAWS 

The curves in fig. [3] suggest the existence of scaling 
laws. In fig. Owe show a plot of the quantity £ = icsL^ 
as a function of the scaled nonlinearity defined by 

9 ^ gL-'/' (9) 

using the same data of fig. [3] One observes, at least for 
the ^ > 1, a clear grouping of the curves, which indi- 
cates that the scaling applies essentially to the nonlinear 
part of the behavior, as intuitively expected: Indeed, as 
4ff(5 = 0) ^ 4, ({g = 0) ~ 96iy-2L-3/'' which is 
dependent of Lq] this fact simply means that Anderson 
localization is not controlled by g. On the nonlinearity- 
dominated region, however, one would (naively) expect 
the nonlinear effects to scale as (w^^) ~ 5-^o ^ for an ini- 
tial state uniformly populating Lq states. Presently, we 
have no convincing explanation for the additional L^^^ 
factor. We note that the scaling is not perfect for the 
small values of Lq, for which the self-trapping is effective 
even for low values of g. 

Fig. [5] shows that the scaling allows us to define 
crossovers between the three dynamical regimes which do 
not depend on Lq. This is confirmed in fig. |51 in which 
we plotted in false-colors the scaled effective wavepacket 
length £ in the parameter plane g, W for Lq = 21 [plot 
(a)] and for Lq = 41 [plot (b)[. Despite of the factor 
2 in the width of the initial state, the two plots are al- 
most identical. The crossover between the quasi-localized 
regime and the chaotic regime is found to be around 
gc ~ 0.1 and the crossover between the chaotic regime 
and the self-trapping regime around g^t ~ 5. 

The use of these scaled variables thus allows us to char- 
acterize the nonlinear dynamics independently of the size 
of the initial state, which constitutes an important step 
in the understanding of these complex dynamics. 
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Figure 8: Survival probability a,t t = IC for states with co- 
herent (empty symbols) and incoherent (full symbols) initial 
phases, for W — 2, Lq = 21 (blue squares) and W = 4, Lq = 
21, (green triangles). The coherent case presents a marked 
enhancement on the transport for g ~ 100. 

V. OTHER EFFECTS 

In this section, we consider three potentially significant 
effects not taken into account above. 

(i) The shape of the initial state. In Fig. [71 we com- 
pared the results discussed above, obtained with the 
"square" initial state defined by eq. ([7|) with results ob- 
tained with an initial gaussian state 

|c„(0)p = ^e-'V2'^^ (10) 

We set (7 = y''(Zo~~17(X(7TT)7l2 so that, for a given 
Lq, both square and gaussian initial states have the same 
second momentum ^ n^|c„p. We see that the dynamics 
is independent of the shape of the initial state to a very 
good approximation. In the self-trapping regime, we see 
that the dynamics depends more on the size Lq than on 
the shape of the initial state. 

(ii) The effect of the quantum phases of the initial state. 
Our choice of using initial random (incoherent) phases in 
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Figure 9: Effect of the sign of the nonlinearity with Lo — 21 
and W — 3. The curves display the survival probability at 
t = 10^ for g > (blue squares) and g < (green triangles). 



ics, which will, most probably, constitute one of the ma- 
jor subjects of atomic physics in the next years. 
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Appendix A: Determination of p{g = 0,t) 



the initial state proved useful in allowing us to give a 
global characterization of the dynamics. However, a par- 
ticular coherent combination of c„ phases can generate 
effects of quantum interference with an impact on the 
dynamics. This is illustrated in Fig. [51 where we com- 
pared the dynamics of an initial state of random phases 
and an initial state with a constant phase for all c„. The 
latter one presents a different behavior, with a marked 
increasing in the transport for g sa 100, which means 
that chaotic behavior is highly favored in this case. Our 
previous result remains however a valid description of the 
"average" dynamics. 

(iii) The effect of the sign of the nonlinearity. In all the 
results presented above g > 0, that is, we have repulsive 
interactions. In the limits of our approach, we have seen 
no significant difference between attractive and repulsive 
interactions, see fig. [HI in contrast with some theoreti- 
cal speculations. This confirms and generalizes a result 
obtained in ref. [4l| in a different context. 

Despite the simplifications we made in order to make 
this problem tractable, it appears that our conclusions 
represent very well the general behavior of the system, 
independently of most microscopic parameters. We can 
thus say that we have characterized the dynamics in a 
rather universal way. 



VI. CONCLUSION 

Nonlinear dynamics is, in general, highly sensitive to 
initial conditions, and no global characterization of the 
dynamics can be made unless one can correctly take into 
account this dependence. This is particularly true for 
the very important system studied in the present work: 
The mean-field generalization of the Anderson model, de- 
scribed by the Gross-Pitaevskii equation. We proposed 
a quantity providing a global characterization of the dy- 
namics independent of the size of the initial state, the 
scaled effective wavepacket length I. This is a first step 
in the necessary development of a coherent language de- 
scribing the effects of nonlinearities in Quantum Mechan- 



The aim of this appendix is to find an analytical for- 
mula describing the evolution of the survival probability 
in the linear case g = 0. Let us first recall a few important 
results in the case of an infinite lattice (no losses - stan- 
dard ID Anderson model). In this case, the correspond- 
ing Anderson Hamiltonian Ha admits eigenstates of 
energy E^. In real space representation, ipi^iji) — {n liy) 
is centered at position Ui, in the lattice, with a localiza- 
tion length £i, < £o{W). As such states are exponentially 
localized in the average, the density of presence of these 
eigenstates at the site n is 

|(n|i/)|2 «tanh(l/4)e-2|"-»-l/^-. (Al) 

In our problem we set an exponential imaginary po- 
tential (absorber) on the border sites of the box Vi (n) = 
—irjaexp {\n\/nc) {\n\ > L/2). In order to make the 
problem tractable, we make the assumption that the 
imaginary potential can be characterize by a constant 
value Vi = ^irja- The corresponding Schrodinger equa- 
tion is 

Decomposing the wave function in the eigenfunctions ba- 
sis, \tp{t)) = '^^aye^^^"*\v), we obtain (setting a„ = 
dajdt): 

ictfi = ^ n) {n\ v) Vi{n)a^ exp [i [E^j, - Ey) t] . 

In the average, the contribution of /i 7^ is negligible 
and, using eq. (jAip . we get 

iau = ^tanh (l/i^) Vi{n) exp (—2 \n — nu\ /iy)^ a^. 

Without loss of generality, we assume rij^ > and note 
£a the distance between the center of the wave packet 
corresponding to this eigenstate and the closest border 
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Figure 10: The blue circles correspond to the values of t,, 
obtained by fitting the numerically calculated evolution of ai, 
with an exponential [eq. (|A2|) ]. The red solid line is a fit of 
these points by eq. (|A3p . 



Figure 11: Graphical determination of p in the weak-disorder 
limit. The curve represents to siv? q and its intersection with 
l{t) determines Qm- The dashed region correspond to the 
integral eq. (|A4[| . 



of the box &t n — L/2. Then 



-rfc tanh(l/£,) ^ e-^^"""")/^" a. 



l + e-2/C^ 



For moderate values of the disorder we can assume 1 4- 
g-2/<'„ ^ 2^ and obtain, integrating the above equation, 



a.(0)e-*/2*'' 
1 



(A2) 



We thus estabhshed two important results: (i) the at- 
tenuation (due to the absorber) of an eigenstate is ex- 
ponential, and (ii) the typical time scale of absorption 
is proportional to exp(2£a/£i^) [l^l- These results can 
be generalized for an imaginary potential that varies ex- 
ponentially. We calculated numerically the evolution of 
Uu (with a non-constant imaginary potential) for a few 
eigenstates (averaged over the disorder) and fitted this 
evolution with equation (|A2p to obtain the correspond- 
ing values of t^, represented as circles in Fig. IIOI plotted 
as a function of in semi- log scale. The dependence 
of ty on can in turn be fitted with the relation (red 
line in Fig. [TO]) 




Figure 12: Comparison between eq. (|A5|) and numerical sim- 
ulations for Lo = 31. The asymptotic behavior of the survival 
probability is very well described by our formula, for different 
values of the disorder, W = 1 (blue squares), W = 2 (green 
triangles), W — 3 (red diamonds). 



with ti, is given by eq. (|A3p . We now take ti, as the 
typical time for decay of the eigenstate v, i.e. we approx- 
imate exp(— t/ty) w 0(ty — t), where 9 is the Heaviside 
function: 

p(5 = o,t) = -l^e(t,-t). 

That is, we consider that the eigenstate is completely 
absorbed at time t^, at which it "touches" the absorbing 
potential. From eq. (jA3|) we see that 



= t2exp[2£a{l/e,-l/eo{W))] 



(A3) 



2e„ 



The fit parameter ^2 thus corresponds to the typical 
time-scale for the attenuation of the eigenstate of largest 
width. 

The evolution of the survival probability p can be 
calculated for a uniform initial state over Lq Anderson 
states, i.e. |a,y(0)p = I/Lq. From eq. (|A2]) 

Pig = 0, i) = ^ |a,|^ = ^ cxp(-tA,) 



log(t.A2)+2^a/^0 

At time t, the surviving states are thus those such that 

^"^^^^"^ log(fA2)+2V^o' 
In the weak-disorder limit one can use eq. ^ and write 

P(5 = 0,i) = / pdq, (A4) 
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where p = I/tt is the density of states [eq. @]. The con- 
dition £q < £{t) sets a maximum value of the quasi- 
momentum that can be determined graphically, as shown 
in Fig. (dll)- The condition 



gives 



and, finally 



^(0 = ^0 sin^ q„ 



- sin-1 m/eo)'^^ 



n((7 = 0, t) = = — sm — ; — 



1/2 

(A5) 



Fitting the numerical survival probability with the 
above expression, best results are obtained by setting 



1 /i _ 4 

2 V 2 4 5 



(A6) 



as shown in Fig. [12] One can think that the distance 
L/2 — io/4 is more likely to represent the typical dis- 
tance from the initial state to the absorber. The factor 
1/2 probably compensates for the brutal step we made 
above of approximating the exponential function by the 
Heaviside function. The presence of the term in ha-s 
a physical meaning: In our simulations, we initially ex- 
cite Lo Wannier states, not Lq Anderson states and as a 
consequence, the initial size of the wave packet is larger 
than Lq. 
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